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Abstract 

CN . The use of Raman scattering techniques to study the mechanical properties of graphene fihxis is 

reviewed here. The determination of Griineisen parameters of suspended graphene sheets under 
"^ ' uni- and bi- axial strain is discussed and the values are compared to theoretical predictions. The 

effects of the graphene-substrate interaction on strain and to the temperature evolution of the 
graphene Raman spectra are discussed. Finally, the relation between mechanical and thermal 
properties is presented along with the characterization of thermal properties of graphene with 
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I. INTRODUCTION 

The growing interest in understanding the mechanical properties of graphene films is 
sparked by the ability to control such properties, and thus to modify the structure and elec- 
tronic behavior for graphene-based applications. Raman spectroscopy is increasingly used to 
measure accurately and non-destructively graphene mechanical or thermal properties, such 
as strain or thermal conductivity. This review outlines the current state-of-the-art in the use 
of Raman spectroscopy to characterize the strain and temperature effects in exfoliated and 
epitaxial graphene. The relationship between strain and film morphology is also reviewed. 

In section II we review the basic atomic structure of graphene, with a brief overview of the 
methods used to isolate and prepare graphene films on various substrates. An overview of the 
mechanical properties of graphene films determined by nanoindentation methods is presented 
in section III, along with the current limitations of such approach. The Raman spectrum 
of graphene in conjunction with its phonon spectrum is described in section IV. A detailed 
overview of the use of Raman spectroscopy for the determination of mechanical properties 
of graphene is presented in section V, with particular emphasis on the characterization of 
strain and of the temperature effects in the graphene films. 

II. GRAPHENE ATOMIC STRUCTURE 

Graphene is a two dimensional sheet of carbon atoms arranged in a honeycomb atomic 
configuration. A single graphene sheet can be folded and, multiple layers can be folded or 
stacked to form sp^ carbon in OD (fuUerenes) , ID (carbon nanotubes, CNT) or 3D (graphite). 
The standard in-plane unit cell of basis vectors lad = jbd = 2.4589 ± 0.0005 A at 297 K^ 
contains two carbon atoms (Fig. la). The resulting two dimensional carbon density is 3.820 
atom- A -. Due to the hybridization of carbon bonds into a sp^ configuration, each carbon 
is bonded to three neighboring atoms in a planar configuration. Two sublattices can be 
identified within a graphene lattice, depending on the orientation of the carbon bonds relative 
to that of their nearest neighbors (Fig. la). The partially filled vr orbitals, perpendicular 
to the graphene plane, are responsible for the electron conduction and the weak interaction 
between a graphene layer and the underlying substrate. This weak interaction is of the 
van der Waals type, independent of the substrate'^ Three possible stacking configurations 





FIG. 1. a) The real space unit cell of a bilayer graphene film in Bernal stacking is shown in gray 
shade, with basis vectors ac and be- Large circles represent the atoms in the top layer, while the 
smaller ones in the second bottom layer. Atoms in grey and black represent the two inequivalent 
sublattices a and /? respectively, b) The reciprocal space unit cell of a single layer graphene, 

— ** 

highlighting the high symmetry points and the reciprocal space unit vectors, a^ and b^j. The first 
Brilloin zone is represented by gray shade. 

exist to form graphitic materials, depending on the relative orientation of the graphene layer 
stacks: Bernal or AB... stacking, hexagonal or AA... stacking, and rhombohedral or ABC... 
stacking. In this review, we only consider the Bernal stacking for multilayer graphene 
films, since the mechanical properties of multilayer graphene have been investigated only 
for this configuration (a comprehensive description of the other stacking sequences can be 
found in ref~). The Bernal (or AB...) configuration is the most common in single crystal 
graphite (80%) by virtue of the lowest stacking energy^. The Bernal configuration is formed 
by stacking two graphene sheets rotated by 60° relative to each other about the z axis. 
The three-dimensional unit cell has 4 atoms, and a third basis vector perpendicular to the 
graphene layer stacks |cg| = 6.672 A at 4.2 K and 6.708 A at 297 K^. The interlayer distance 
is cg/2. Because of the 60° rotation between the subsequent layers, the two sublattices 
in graphene see a different local environment in the Bernal configuration: an a atom is 
positioned directly above an a atom in the sheet below, whereas a /3 atom is positioned 
above the (empty) center of the ring of the sheet below (Fig. la). The presence of a non- 
graphitic substrate alters the equivalence between the two sublattices with possible effects 
on both the mechanical and electronic properties, as discussed in Section V. 



The Brillouin zone for a single graphene layer is shown in Fig. lb. It exhibits high 



symmetry points: the F point at the zone center, the M point in the middle of the hexagonal 
sides and the K and K' points at the corners of the hexagons. K and K' are inequivalent 
points, since they correspond to the two different and inequivalent sublattices in the graphene 
atomic structure. 

Graphene samples can be prepared by mechanical exfoliation of highly oriented pyrolithic 
graphite (HOPG) '"' , which leads to the production of micrometer scale single and multilayer 
graphene sheets with high degree of control over their thickness. Graphene can be also 
grown epitaxially on SiC surfaces by high temperature Si sublimation, in ultrahigh vacuum 
(UHV)'"*'" and in controlled environment ^"^'"^^ Epitaxial graphene can also be grown on the 
surfaces of various metals such as Pt^\ Ni^"'^'\ Ir^''^'*^, Ru^*^*"^^ and Cu"". With this method, 
large domains can be obtained (domain size ~ lO/im)^ '. Epitaxial graphene grown on metals 
can be transferred from the synthesis substrate to any chosen substrate' '. This procedure 
is suitable for investigation of large scale graphene layers either suspended or transferred 
to various substrates. The graphene-substrate interaction strongly depends on the type of 
substrate due to the different degree of adhesion of graphene to the substrate (whether, 
for example, graphene is grown epitaxially on a substrate or mechanically transferred to 
it). Therefore the choice of substrate and synthesis method have several implications in the 
mechanical properties of the epitaxial graphene film. 

III. GRAPHENE MECHANICAL PROPERTIES MEASURED BY NANOINDEN- 
TATION 

The mechanical behavior of graphene layers can be described macroscopically by contin- 
uum elasticity theory. In this spirit, nanoindentation techniques are well suited to measure 
the macroscopic mechanical properties of graphene, including Young's modulus and bending 
stiffness. For example, by using nanoindentation methods on suspended multilayer graphene 
flakes, the bending stiffness has been measured and found to be in the range from 2 x 10~^^ 
N/m to 2 X 10"^^ N/m for 8 to 100 layers, respectively. Static nanoindentation experiments 
based on the deflection of AFM cantilevers pressed within 100 nm of the center of ~ 1/im 
long double-clamped graphene fllms, provided a measurement of the effective spring con- 
stant of multilayer graphene (1-5 N/m). The spring constant was found to scale with the 
dimensions of the suspended region and the layer thickness (from 5 to 30 layers), and of the 



extracted Young modulus of 0.5 TPa, independent of thickness^^ A significant limitation of 
the use of nanoindentation techniques is the requirement of a graphene layer to be suspended. 
The presence of a substrate, over which graphene may either be deposited (Si02'', glass and 
sapphire"'' or polymers^"''""') or directly grown epitaxially (e.g. SiC^'' and metals^"""^), makes 
it hard to separate by nanoindentation measurements the intrinsic mechanical properties of 
a graphene from that of the substrate. 

In contrast to nanoindentation, Raman spectroscopy provides access to information relat- 
ing to the underlying chemical bonds. Besides complementing the coarse grained approach of 
macroscopic elasticity, the interrogation of bond vibrations by optical spectroscopy enables 
the retrieval of information about mechanical and structural properties of films that can have 
monolayer thickness and be strongly interacting with a substrate. Raman spectroscopy has 
thus been used to measure mechanical properties of graphene films, both freestanding and 
on a substrate^'"'''^'', at room and at elevated temperatures^^''^'^. 

IV. RAMAN SCATTERING IN GRAPHENE AND GRAPHITE 
A. Raman Spectroscopy of Graphene 

The Raman spectrum of carbon based materials is characterized by a set of common 
features in the region between 800 and 2000 cm~^, in particular the so-called D and G bands, 
which lie at around 1330-1360 and 1580 cm~^ respectively for visible excitation^'^''^^'^'^ as 
shown in Fig. 2a. Under these excitation conditions, the Raman spectra of carbon films are 
dominated by the sp^ sites, because visible excitation always resonates with the vr states. 
Due to the comparatively small cross-section for the amorphous sp^ vs sp^ C-C vibrations, 
a significant fraction of sp'^ bonds is required in a sample for the sp^ peak at 1332 cm~^ to 
be visible, as is the case in diamond (Fig. 2)-". 

The phonon dispersion curves of graphene (Fig. 3) are the key to understand its Raman 
spectrum. They consist of three acoustic phonon modes (A) and three optical (O) phonon 
modes since the graphene unit cell contains two carbon atoms (Fig. la). Among these modes, 
one acoustic branch and one optical phonon branch correspond to out-of-plane phonon modes 
(o), while for the other acoustic and optical phonon branches, the vibrations, and thus the 
phonon modes, are in-plane (i). Each in-plane mode has two branches, one longitudinal (L) 
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FIG. 2. a) Raman spectra of graphite^^, single-layer graphene'^" , metallic and semiconducting 
carbon nanotubes^'', low and high sp'^ amorphous carbons^'', and diamond'' for visible excitation 
(excitation energy: 514 cm^-*^). b) "Molecular pictures" of the the Fi2g and Ai^ modes, correspond- 
ing to the G and D peaks, respectively. 

and one transverse (T). Following the high symmetry FM and FK directions, the six phonon 
dispersion curves are assigned to LO, iTO, oTO, LA, iTA, and oTA phonon modes '"''"'. In 
graphite the LO and iTO modes are degenerate at the center of the Brilloin zone, the F point. 
According to group theory, these modes are the only Raman active modes, corresponding 
to the two dimensional E2g phonon. The G peak (located around 1580 cm~^) corresponds 
to such doubly degenerate E2g mode at the Brillouin zone center'^^''^^ In the "molecular" 
picture of carbon materials, the G peak is due to the bond stretching of all pairs of sp^ 
atoms (Fig. 2b). 

The D peak (~ 1340 cm~^) corresponds to modes associated with transverse optical 
(iTO) phonons around the edge of the Brillouin zone (K or Dirac point)-". In the molecular 
picture, it is associated with the breathing mode of the sp^ aromatic rings (Fig. 2b)'^^''^^. 
The D peak is energy dispersive, so that its position is dependent on the excitation energy 
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FIG. 3. Phonon dispersion plot of a single-layer graphene, calculated (lines)"' and experimental 
(points)'"'. Different experimental points corresponds to the different branches. 



(Fig. 4)'^''*. The D peak is usually very intense in amorphous carbon samples, while it is 
absent in perfect graphitic samples. Its overtone (2D, ~ 2660 — 2710 cm~^) however is 
always visible even when the D peak is absent. Such peculiar behavior is due to the double 
resonance (DR) activation mechanism ■'' of the D peak, which requires the presence of defects 
for its initiation'"'"'. In a double resonance process, Raman scattering is a four-step process: 
(i) a laser induced generation of an electron-hole pair; (ii) electron-phonon scattering with an 
exchanged momentum q~K; (iii) electron scattering from a defect, whose recoil absorbs the 
momentum of the electron- hole pair; (iv) electron- hole recombination-'". The requirements 
of conservation of energy and momentum can only be satisfied if a defect is present. In a 
perfect sample, momentum conservation would be violated by the DR mechanism, and thus 
the D peak is absent. Momentum conservation however is always satisfied in case of the 2D 
peak, without the need for defect activation, since the process involves two phonons with 
opposite momentum vectors^"'. A similar process is possible with scattering within the same 
K point. This intra- valley process activates phonons with small momentum q, resulting in 
the so-called D' peak, located around ~ 1620 cm~^ in defective graphite ^^. 

Scattering from holes can also occur in the Raman process. In graphene, under these 
circumstances, the electron is not scattered back by a phonon of momentum -q, but instead 
a hole is scattered forward by a phonon with momentum +q. In this case, during the 



electron-hole generation, both electron and hole scattering processes are resonant. The 
electron-hole resonant recombination at the opposite side with respect to the K point is also 
resonant, resulting in the triple resonance scattering process (TR). It has been suggested 
that the higher intensity of the 2D peak relative to the G band in a graphene monolayer is 
due to the triple resonance activation mechanism'^'^ 

B. Graphene Metrology with Raman Scattering 

Raman spectroscopy, as a non-invasive probing technique, has been extensively employed 
to characterize graphene layer thickness'^*^''^, domain grain size^'^''^^^'^'^ doping levels'^'^'**"'^^, 
the structure of graphene layer edges^'^""^, anharmonic processes and thermal conductivity^^'^'^. 
This has been possible through a combined investigation of the Raman peaks D, G and 
2D in graphite and graphene films of various thicknesses and morphologies. An indicative 
comparison of the Raman spectra of graphene and bulk graphite is made in Fig. 4a'"\ The 
most striking difference between the individual graphene layers and graphite resides in the 
change in shape and intensity of the 2D peak. While the 2D peak in graphite consists of two 
peaks 2Di and 2D2 (with intensities of 1:4 and 1:2 compared to the G peak, respectively), 
the 2D peak in one single graphene layer has only one component with roughly four times 
the intensity of the G peak (Fig. 4a). For multilayer graphene (Fig. 4b), the evolution in the 
shape of the 2D peak has been used to determine the layer thickness'^*^*'^^'^'*. The splitting 
of electronic bands in bilayer graphene is responsible for the splitting of the 2D peak into 
four components"' (Fig. 4c). The two lower components further decrease while the higher 
wavenumber components increase as the film thickness approaches 5 layers. Above this 
threshold, however, the determination of the layer thickness with Raman becomes rather 
difficult, as the shape of the 2D peak is increasingly similar to that of bulk graphite. 

Early investigations of disorder in graphitic carbon'^^ show that the ratio of the D and G 
band intensities (Id/^g) is inversely proportional to the in-plane crystallite size Lq, measured 
independently with x-ray diffraction. Such relation, known as the Tuinstra-Koenig (TK) 
relation, has been refined in recent years to provide an empirical method to determine the 
size of graphene domains from the Raman spectrum under a given excitation energy'^^''^'^. 
There are known limitations in this approach, as the distribution of domains with different 
sizes is such that the smaller domains are weighted more, leading to an underestimation of 
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FIG. 4. a) Raman spectra of bulk graphite and single layer exfoliated graphene, taken with 
excitation energy 514 nm. b) Evolution of the Raman 2D spectra with layer thickness, taken with 
laser excitation 514 nm. c) The 2D peak in the graphene bilayer is composed of four Lorentzian 
components, while the single layer has only one. The dispersive nature of the 2D peak is clearly 
visible in the net shift of the 2D peak in plots c, with excitation A = 514 nm, when compared to 
its position for an excitation of A = 633 nm''". (Reprinted with permission from ref."*". Copyright 
2006 American Physical Society.) 



the average size distribution. In addition, the use of peak intensity ratio instead of peak area 
ratio, underestimates the average domain size, since the full-width-half-maximum (FWHM) 
of the D peak increases significantly in comparison to that of the G peak-". Furthermore, 
the ratio Id/Ig is known to depend on the electron concentration (and thus on the film 
doping)^'', limiting the application of the TK relation when the doping concentration is 
unknown. Regardless of the limitations, the use of the TK relation allows an estimation of 
the degree of disorder in the graphene film. 



V. PROBING MECHANICAL PROPERTIES OF GRAPHENE WITH RAMAN 
SPECTROSCOPY 

Any changes in the atomic structure in a crystalhne sohd due to plastic deformation, 
strain, or thermal expansion are reflected in the phonon spectrum of the crystal. By prob- 
ing the phonon spectrum with Raman spectroscopy, such changes can be detected, thus 
providing insight into the mechanical and thermal properties of materials such as graphene. 

Strained semiconductors have received signiflcant interest in the past because of the wide 
ranging implications of strain, such as the ability to engineer the electronic structure and to 
affect the carrier mobility in silicon-based materials for electronic device application ■'^"^. The 
application of an external stress on a crystal results in a lattice strain, i.e., in a change in 
interatomic distances and consequent redistribution of electronic charge. Isotropic compres- 
sion (hydrostatic pressure) generally results in an increase in the frequency of the vibrational 
mode (phonon hardening), while isotropic tension results in the decrease in the vibrational 
frequency (phonon softening). Application of anisotropic stress has more complex effects, 
and can result in lifting of the degeneracy of phonon frequencies. 

In graphene, changes in the Raman spectra have been observed as a consequence of 
the presence of stress, either induced artiflcially on suspended or exfoliated graphene^"' ■'''^"'''^ 
or provided by the interaction with the substrate for graphene grown epitaxially on SiC 
substrates^'^'^'^'^^"^'^. Such changes consisted of a systematic upshift in the position of the 
main Raman D (when present), G, and 2D peaks, by up to 30, 31 and 64 cm~^, respectively^^, 
for an applied strain of up to 1.3 %. 

When a uniaxial tensile stress is applied to a graphene layer, the splitting of the G peak 
has also been observed, reaching up to 15 cm~^, for an applied strain of 1.3%"^ '''''. Each peak 
in the split G band corresponds to two orthogonal modes, having eigenvectors perpendicular 
to the applied strain (E^ ) and parallel to it (E^ ). When the uniaxial compressive strain 
is applied, sp^ bonds along the direction parallel to the applied strain are shortened and 
hardened, while those perpendicular to it are only slightly affected (Fig. 5). Hence, under 
uniaxial strain, only the peak G~ corresponding to the E^ mode is signiflcantly shifted 
relative to the unstrained E^g (by as much as 30 cm~^), while the peak G^ corresponding 
to the E^ mode is only moderately shifted (up to 15 cm^^). Since this effect is purely 
mechanicaP^, the full- width-half- maxima of G~ and G^ remain constant. The FWHM of 
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FIG. 5. Raman spectra in the single-layer graphene of the a) G and b) 2D peaks, as function of 
the applied uniaxial compressive strain percentage, indicated in the right side of each spectrum. 
The spectra are acquired "with the polarization of the incident light parallel to the direction of the 
strain. The double degeneracy of the G band is broken as a consequence of the applied strain, 
resulting in t"wo peaks G"^ and G~"''. (Reprinted with permission from ref.-'. Copyright 2009 
American Physical Society.) 



the 2D band is also unchanged. A similar behavior is observed in carbon nanotubes, where 
the tube curvature induces the splitting of the G band, with a significantly larger shift for 
the component parallel to the curvature'". The intensities of the two peaks G~ and G"^ vary 
with the polarization of the scattered light along the direction of the strain, allowing the 
sample crystallographic orientation with respect to the strain to be probed^^'^^. 



In spite of specific changes in the electronic and vibrational band structure, the strain- 
induced frequency shifts of the Raman active E2g and 2D modes are independent of the 
direction of strain, which has been observed experimentally"'^ and confirmed by ab initio 
calculations'^^. Thus, the amount of strain can be directly determined from a single Raman 
measurement^^ . 
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A. The Griineisen Parameter for Uni- and Biaxial Strain 

The rate of change with strain of a given phonon frequency in a crystal is determined 
by its Griineisen parameter'^ '''^*^. In metrology applications, accurate values of Griineisen 
parameters are crucial for quantifying the amount of strain in the system, reflected in the 
change in phonon frequency from its value in the absence of strain. In presence of uniaxial 
strain, the Griineisen parameter for a particular band m associated with in-plane Raman 
active phonon band (where m is either the D or G band in graphene), is defined as*^': 

where et = ^n + etj is the hydrostatic component of the applied strain with / and t referring 
to the directions parallel and perpendicular to the applied strain respectively, and cu^ and 
u^^ correspond to the phonon frequencies of peak m at zero strain and in presence of an 
applied strain, respectively. For a given shear component of strain, e^ = en — eu, the shear 
deformation potential f3m is defined as: 

f^m = V-^ (2) 

For the G band corresponding to the E2g phonon, the shifts in the two components G"*" 
and G~ relative to the position at zero strain, ujq, are given by: 

where Auq and Auq are the shifts associated with the hydrostatic and shear components 
of the strain respectively. Under condition of uniaxial strain, eu = e and e^ = —ue, where 
z/ is the Poisson ratio'" . In case of graphene, if the layer adheres well to the substrate used 
for strain analysis, such as for example polyethyleneterephtalate (PET)-'^, the Poisson ratio 
of the substrate must be used, instead of in-plane Poisson ratio for bulk graphite. Under 
uniaxial strain. Equations (3) can be solved, yielding both the Griineisen parameter and the 
shear deformation potential for the G band, as functions of the shifts in the positions of the 
two components G"*" and G~: 
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po = - ^:°r/"° (5) 

Under the conditions of biaxial strain, en = e^ = e, there is no shear deformation potential 
and no splitting of the G peak. In this case, Equation 3 can be solved to provide the 
Gruneisen parameter^"''-*^: 

It is however possible that local anisotropics in the applied biaxial strain, possibly induced 
by the substrate over small domain size (such as in epitaxial graphene grown on SiC), may 
cause an increase in the FWHM as a result of a local splitting of the G band. It is also worth 
noting that under biaxial strain conditions, the shift in the peak position is independent of 
the presence of any substrate, because of the absence of a sheer deformation term and thus 
the absence of the Poisson term u in eq. (6)"'. 

The Griineisen parameter can be similarly derived for the D and D' bands in graphene. 
Of the two, only the first is single-degenerate, and corresponds to Ai^ phonons at the K 
point (Fig. 3). The D peak is thus not expected to split under uniaxial strain, and only the 
hydrostatic component of the stress is present. The Griineisen parameter for the D peak 
(which is equivalent to that of the overtone 2D) can be written as: 

^uniax _ AcJd,2D /^x 

or: 

(Note that the shear deformation potential (3 for the D and D' bands cannot be extracted, 
because of the lack of shear component of the applied uniaxial strain). The D' band is 
associated with an E symmetry mode, which is double-degenerate; as such, a splitting is 
expected under uniaxial strain. Experimentally the only study to report on the effects of 
strain on the D' peak did not observe any splitting, due to the weak intensity of this peak 
and the small range of applied strain" '. For small strains, the Griineisen parameter for the 
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TABLE I. Griineisen parameter and shear deformation potential for a single layer graphene. The 
Griineisen parameter 72D is equivalent to that of 7^). * 721) has been measured directly only in 
case of biaxial strain ' - . 

D' follows eq. (8). In the case of biaxial strain, equation (7) is the same as equation (4), 
which can be generalized as: 



hiax 
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Aa;„ 



(9) 



where m corresponds to the D, G or 2D bands. It is worth mentioning that in all cases, 
the detection of strain effects is the most sensitive if the 2D band is considered. With a 
spectrometer resolution of ~ 2 cm~^, the sensitivity for uni- and biaxial strain is 0.03 and 
~ 0.01, respectively. 



B. Determination of the Griineisen Parameter in Graphene 



Mohiuddin et al. provided a complete characterization of the Griineisen parameters for 
the G and 2D bands of exfoliated graphene^ '. In order to measure the Griineisen parameters 
and the shear deformation potential of a single layer exfoliated graphene, equations (4), (5) 
and (7) were used to fit measured shifts in the positions of the G and 2D Raman bands 
as a function of applied uniaxial stress. The resulting Griineisen parameters and shear 
deformation potential for the main Raman peaks for a single layer graphene are summarized 
in table I, along with previous theoretical and experimental studies. The values of 7^ and 
Pg from ref.""' are in good agreement with those calculated with density-functional theory 
(~ 1.8''^) and first principle calculations (1.87"'). Metzger et al. measured directly the 
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Griineisen parameters for biaxial strain, by placing a single-graphene layer onto a substrated 
prepatterned with shallow depressions*^". The adhesion of graphene to the substrate across 
the pre-patterned depression, despite the induced biaxial strain, allowed a controlled and 
precise determination of the biaxial strain. By using eq. (9), the Griineisen parameters for 
the G (2.4) and D (3.8) peak were extracted and found to be higher than those measured from 
uniaxial strain-' or calculated""''*'^. It has been speculated that the larger values of both the 
Griineisen parameters and the peak shifts when compared to previous measurements were 
due to a better adhesion of graphene to the substrate in the latter studies. This leads to 
a measurement of the strain actually transferred to the graphene layer from the substrate 
(i.e., with no slippage). However, the large difference between the measured values of the 
Gruneisen parameters for the G band (1.8-"' vs. 2.4*'") but not for that of the D peak (3.55^^ 
vs. 3.8*'*') remains unexplained. Recently, Tsoukleri et al. investigated the effects of slippage 
of the graphene layer from the substrate by measuring the strain in a graphene layer, either 
supported on a poly(methyl methacrylate) (PMMA) cantilever beam or embedded into the 
cantilever' '^ The values of diJ^""^ j de were measured by altering the stress applied to the 
cantilever. The measured values for either the supported or embedded case are in good 
agreement with those measured initially by Mohiuddin et al' ', suggesting a negligible role 
of slippage. When compressive stress is applied, however, the effects of buckling and the 
formation of ridges in the supported graphene layer determine a significantly different value 
of duj2£}"'^ /de, compared to the embedded case. 

By applying such parameters to eq. (4), the gradients in Raman peak position per unit 
of applied strain are extracted. A summary of both theoretical and experimental studies is 
reported in table II. The use of the correct value for the Griineisen parameter is extremely 
important, because it affects the estimated value of the Raman peak shift for a given strain. 
Often the Griineisen parameter of CNT is used, leading to a questionable estimate for the 
gradient in the peak position. For example, 72D = 1.24 '*' duj^""^ j de ~ —27.1 cm~^/% '*' to 
be contrasted to ~ —83 cm~^/% when ^^d = 3.55 is used per ref.-"', obtained on a single- 
layer graphene. This result has been used to justify the measured value of the gradient 
in peak position for uniaxial strain. However, the absence of any splitting of the G peak 
and lack of any difference in Raman peak position between uni- and biaxial graphene"'*', 
which are in contradiction with theory, suggests that the applied strain is either far from 
being uniaxial"'*' or points to poor sample quality. As a further indication, the estimated 
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TABLE II. Gradients in Raman peaks position per units of applied strain (cm~^/%), for a single 
layer graphene. Whenever the D peak was not present in the measured Raman spectra of the single 
layer graphene, the gradient in the shift of the D peak is taken as half that of its overtone, the 2D 
peak, as expected from eq. (7) and (9). * Measurements carried out under compressive strain. 

gradient in the G peak position as function of applied strain [duc/de ~ —27.8 cm~^/%) is 
consistent with the averaged value of the gradients of the shifts in the G^ and G~ peaks 



[Oooq /ae 



-27 cm 



'"). This is also consistent with the average value obtained 



from measurements on carbon fibers (~ —25 cm~^/%), where individual sub-bands cannot 
be distinguished due to the broad G band for amorphous carbon'''. The similarity in such 
measurements between graphene and graphite indicates that the in-plane Young modulus 
for graphene and bulk graphite are similar' "\ 

C. Substrate-Induced Strain on Graphene 



While uni- and bi-axial strain can be artificially applied to suspended graphene layers, 
strain can arise in graphene heterostructures from the interaction between graphene layers 
and the underlying substrate. Initially, in the case of exfoliated graphene, no appreciable 
shifts were observed in the G band of a graphene layer transferred onto Si02/Si and GaAs 
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Distance Oim) 

FIG. 6. a) Atomic force micrograph of a graphene layer suspended over a microfabricated trench, 
after thermal cycling to 700 K. b) Spatial mapping of the G band Raman shift taken perpendicular 
to the trench before and after the thermal cycling to 700 K '''. (Reprinted with permission from 
ref.'"'*^'. Copyright 2009 American Chemical Society.) 

substrates^ ^ Small downshifts of about 5 cm~^ were observed for graphene placed on sap- 
phire and glass, with a splitting of the G band in the latter case"^ The higher adhesion 
offered by sapphire substrates, which is sufficient to introduce a small amount of strain in 
the graphene layer during the mechanical placement, was attributed to the particular bind- 
ing of carbon-sapphire^^. The binding is also responsible for the growth of highly aligned 
CNT on sapphire substrates'''. Recently, however, more comprehensive investigations of the 
evolution of the mechanical and morphological properties in graphene suspended over a mi- 
crofabricated trench reported variations in the positions of the G and 2D bands, during and 
after thermal cycling"^*^. Upon thermal cycling to 700 K, while in purely suspended regions 
no shifts were observed, large upshifts (~ 23, ~ 10, ~ 5 cm~^ for a single-, bi- and tri-layer 
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FIG. 7. Raman spectrum of the Si-terminated SiC clean surface is compared with that of 1.5 
epitaxial graphene layers. The G peak of graphene (indicated in correspondence to the dotted line, 
at ~1592 cm~^) is convoluted with the second order peaks of the SiC substrate. The scattering 
contribution of the SiC substrate can be removed by using a depolarized scattering configuration 
(as shown with the solid lines) ^''. Excitation energy: 633 nm. (Reprinted with permission from 
ref.'^''. Copyright 2008 American Physical Society.) 

graphene respectively, corresponding to a compressive strain of 0.39%, 0.18% and 0.09% 
respectively) were instead observed in the regions where graphene was in contact with the 
underlying substrate. While graphene was compressed in the region over the substrate, the 
compression was relieved and the formation of ripples was observed in the purely suspended 
region (Fig. 6)''''^. 

The role of the substrate on strain in graphene films has been also investigated ex- 
tensively on graphene grown epitaxially on SiC surfaces (so called epigraphene) by high- 
temperature decomposition^^''^'^'®^"^^ Figure 7 shows representative spectra of a single crystal 
6H-SiC(0001) surface and that with 1.5 layers of epigraphene. The Raman peak of zone- 
center optical (G) phonons in monolayer epigraphene is overwhelmed by the second order 
signal from the SiC substrate, a broad band occupying the same spectral region. This un- 
fortunate coincidence limits the ability to measure precisely the position of the epigraphene 
G band itself. This limitation can be overcome by the use of a depolarized scattering 
configuration^*^, as shown in Fig. 7. 

Raman spectra of epigraphene on the Si-terminated 6H- and 4H-SiC (0001) substrates 
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FIG. 8. Raman spectra of single- and bilayer epitaxial graphene on Si-terminated SiC, SiC sub- 
strate, micromechanically cleaved/exfoliated graphene (MCG) and bulk graphite as indicated. The 
shift in the position of the 2D peak is shown in the inset. Excitation energy: 532 mn '"'. (Reprinted 
with permission from ref. '^. Copyright 2008 American Physical Society.) 

usually show a blueshift in the graphene epilayer peak positions, with respect to those on 
exfoliated graphene''"'''-. The extent of shift is different for each Raman peak, as shown 
in Fig. 8. The blueshift recorded varies by up to 22 cni~^ for the G peak and 64 cm"^ 
for the 2D peak. Graphene on C-terminated SiC substrates have not been investigated in 
full details with Raman spectroscopy. However it is speculated that the decoupling of the 
graphene layer grown on the C-termination may reduce the amount of strain in the film. 

The large shift in epitaxial graphene layers on Si-terminated SiC was attributed to com- 
pressive strain in the graphene layer. This explanation may seem surprising, since no external 
strain was applied to the system. However, the only possible alternative explanation, charge 
transfer from the substrate, was ruled out, based on the fact that it could not account 
for the magnitude of the shifts in the G and 2D peaks. Indeed, while charging induces a 
shift in the G peak up to ~20 cm^^ for an electron concentration of 4x10^'^ cm~^''', the 
shift in the G peak corresponding to charge measured in a monolayer graphene on 6H-SiC 
(1.4x10^^ cm~^") would only account for approximately 7 cm~^. Similarly, shifts in the 2D 
band corresponding to the given amount of charge in monolayer graphene is negligible ^*^. 
Hence the observed shifts could only be explained in terms of strain in the system^^'^^'^^'^^. 
By using the Griineisen parameters evaluated under applied uni- and biaxial strain on sus- 
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FIG. 9. Coefficient of thermal expansion as a function of temperature for a single layer of graphene 
(calculated)*"^, SiC (measured)'^, and graphite. Note that the coefficient of thermal expansion of 
graphene is always negative. 

pended graphene layers (Table I), the amount of intrinsic strain in epigraphene can be 
evaluated using equation (9). It is interesting to note that the shifts of the D and G peaks 
occur in the approximate ratio of 1 : lA-^'^~, which is in good agreement with the ratio be- 
tween the Griineisen parameters for those peaks on exfoliated graphene in presence of biaxial 
stress (1.8 : 2.7, table I). Hence, for the maximum observed upshift of 22 and 64 cm~^ for 



the G and 2D peaks, the corresponding strain in epigraphene is approximately 0.7-0.8^ 



'26 



The shifts in the Raman spectra are found to decrease as the number of graphene layers 
increases. More specifically, the G and 2D peaks in the epitaxial graphene bilayer are found 
to be shifted by up to 7 and 22 cm~^ (as opposed to 22 and 64 cm~^ for the monolayer, 
respectively), to approach the unstrained values for films thicker than a ~6-9 layers'^-^. 

The presence of strain in epigraphene was initially explained in terms of the differ- 
ence between the lattice constant of the reconstructed 13x13 graphene layer supercell 
(ctgraphene = 31.923 A) and of the reconstructed SiC 6v^x6v^ supercell {astc = 31.935 A)''^. 
Such small difference cannot account for the significant amount of strain measured. Com- 
pressive strain at room temperature in the graphene layer was later attributed to the large 
difference in the coefficients of thermal expansion (Fig. 9) between graphene (ogr, as mea- 
sured and calculated in ref.*^*^) and SiC (asic? as measured in ref. '"*) during cooldown from 
the synthesis temperature"'^'''^. This difference Aa(T) is nearly constant between room tem- 
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perature (RT) and the graphene synthesis temperature, Tg ~ 1250°C. If the epitaxial film 
is in mechanical equilibrium with the SiC surface, as a stress-free monolayer commensurate 
with the 6 X -\/3-reconstructed SiC surface at T5, a large compressive strain would develop 
in the film upon cooling, since SiC contracts on cooling, while graphene expands^'': 



1 



1 

= exp 



dT'Aa{T') 

RT 



(10) 



Ferralis et al. found that the shift observed in the position of the Raman peaks strongly 
depends on the duration of the high temperature annealing^^'''^. The evolution in the shift 
of the 2D peak as a function of the annealing time is shown in Fig. 10. It was observed 
that for short annealing times (up to 2 minute) the G and 2D Raman peaks were almost 
unshifted from their unstrained values. Longer annealing times (up to 1 hour) were found to 
produce the largest shifts (as high as 22 cm~^ for the G band, corresponding to a strain of 
~0.8%, based on eq. 9). It was argued that higher compressive stress at room temperature 
resulted from a lower stressed film at the synthesis temperature {Ts)., while a nearly stress 
free film at room temperature indicated that the film existed under high tensile stress at 
Ts- Within experimental accuracy, the strain measured at room temperature might well 
vanish for very short annealing times. In contrast, for long annealing times, the graphene 
layer reaches mechanical equilibrium with the substrate at the synthesis temperature T5, 
and a compressive strain develops at room temperature film (up to ~0.8%). This analysis 
suggests that mechanical equilibrium with the 6--\/3 SiC substrate at Ts is indeed achieved 
for annealing times longer than 10 minutes, while for shorter annealing times (~5 minutes 
or less), graphene is under high tensile strain at T^^^'*^^. 

A direct correlation between the strain distribution and graphene surface morphology 
was made using a combined Raman spectroscopy and electron channeling contrast imag- 
ing (ECCI)*^^. It was found that the roughness of the SiC substrate terraces from where 
epigraphene grows increased paralleling the increase in the Raman peak shifts under the 
same conditions, as shown in Fig. 10. This observation provides a possible mechanism for 
strain relaxation. For long enough annealing times, tensile strain developed at Ts is relieved 
by the roughening of the step edges to which graphene films are pinned. Such increase in 
roughness does not induce a significant change in surface coverage (±0.2 ML). For short 
annealing times, surface relaxation and roughening do not take place, leaving the SiC ter- 
races morphologically unchanged. Similarly, large inhomogeneities in the distribution of 
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FIG. 10. The evolution of the shift in the 2D peak as a function of the anneahng time is com- 
pared with the evolution of the edge roughness Redge of the Si-terminated SiC terraces after 
graphitization''^ . The edge roughness Redge is defined as the difference in the normalized aver- 
age mean square deviation of any graphitized terrace edge with that of the initial ungraphitized 
surface. Several profiles of terrace edges are extracted from electron channeling contrast images of 
samples prepared with the same annealing time and temperature. Each profile (black curve on the 
ECCI image of a sample annealed for 8 minutes) is fit with a 9^^ order polynomial to obtain an 
edge baseline. The normalized average mean square deviation (and thus the edge roughness Redge) 
is extracted from the baseline. (Reprinted with permission from ref.'^^. Copyright 2008 American 
Institute of Physics.) 



strain within the same epigraphene layer were reported by combined Raman mapping and 
atomic force microscopy (AFM)*^^ Large shifts in the 2D Raman band (up to 74 cm~^, 
corresponding to a strain of about 1.0%) were observed to correspond to regions with screw 
dislocations, step terraces and macrodefects, while regions with less pronounced band shifts 
corresponded to large fiat terraces (Fig. 11). The strain distribution map obtained with 
Raman spectroscopy appears to be correlated with the surface morphology of the graphene 
film, monitored by AFM, confirming that changes in the physical topography are related to 
changes in the strain of the graphene film^'^ . 
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FIG. 11. Raman spectral map corresponding to the position of the 2D peak of epitaxial graphene 
grown on the Si- face of SiC. (a) is compared to the AFM image (b) of the same area, (a, b) near a 
SiC screw dislocation where its position is marked with an "x" and (c, d) where such defects are not 
present. The Raman topography is correlated with the surface morphology of the graphene film as 
revealed by AFM, suggesting that changes in the physical topography are related to changes in the 
strain of the graphene film''''. (Reprinted with permission from ref.*'^ Copyright 2009 American 
Chemical Society.) 

D. Characterization of Thermal Properties of Graphene with Raman Spec- 
troscopy 



Raman spectra have a significant temperature dependence, both in intensity and in posi- 
tion of the Raman peaks. For example, the ratio of the intensities of antiStokes and Stokes 
peaks is commonly used as a metrology tool to determine the actual temperature of the an- 
alyzed sample""". Since strain in the lattice also affects Raman peak positions, it is crucial to 
understand and discern the role played by the changes in lattice parameters (due to strain 
or thermal expansion) from purely isovolumetric thermal dependencies. Experimentally, 
separating the two contributions is complicated, especially if either mechanism is not easily 
controllable, or strictly depends on the position of the Raman peak for its determination. 
In complete absence of strain, shifts in Raman peaks observed in response to temperature 
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changes reflect both elementary anharmonic processes (electron-phonon and phonon-phonon 
scattering) and changes in lattice parameters with temperature (thermal expansion). The 
temperature dependence of the G and 2D peak positions Um for single and bilayer suspended 
graphene is approximately described by^^'^®: 

^m = ^m + XmT (11) 

where cu^ is the position of the peak m (either G or 2D) at T=0 K, and Xm is the first- 
order temperature coefficient of the same peak. By measuring the position of the G and 2D 
peaks as a function of sample temperature, the temperature coefficients are extracted for 
single and bilayer graphene (Fig. 12). The results are reported in table 111, and compared 
with other carbon-based materials. It should be noted that the geometrical configuration 
employed in these experiments (a graphene sheet rigidly connected to the substrate) does 
not guarantee the conditions of a strain- free environment. Hence, the actual determination 
of the thermal evolution of the Raman spectrum through these experiments may include 
non-negligible contribution from strain. 

The temperature dependence of the G peak for the single layer is found to be higher than 
for the bilayer. Both values are higher than that for HOPG, and are expected to approach 
the HOPG value for thicker graphene films. The temperature coefficient Xm depends on the 
anharmonic potential constants, the phonon occupation number and the thermal expansion 
of the graphene two-dimensional lattice^^. The contribution of anharmonic terms is most 
significant at high temperatures; hence the overall thermal dependency is not expected to 
follow a linear trend"'^. The non-linearity must be taken into account when using calibration 
of thermally induced shifts in the Raman spectra of graphene. Commonly used linear fits 
need to be accompanied by the temperature range used for the measurements, as reported 
in Table 111. In HOPG, Xm is found to depend mostly on the anharmonic contribution, due 
to direct coupling of phonon modes. Since thermal expansion occurs primarily along the 
c-axis, its effect on the in-plane G and 2D Raman modes are not very pronounced'^''. 

It is however important to note that the interaction with the substrate may strongly affect 
thermal expansion of graphene, resulting in a different value of Xm for purely suspended 
versus strongly interacting graphene layers. This might be the cause of the different value 
measured by Cai et al. for a single layer graphene grown by CVD and pressed against a 
Au/Si02 thin film on Si'"' . As a further indication of a strong interaction with the substrate, 
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FIG. 12. Temperature dependence of the G peak position (shown in the inset) for the single- 
(a) and bilayer (b) exfohated graphene. The measured data is used to extract the temperature 
coefficient for G peak-''^*^. (Reprinted with permission from ref."'. Copyright 2007 American 
Chemical Society.) 

the same value of x ^^"^ found on regions of the same graphene layer either supported and 
suspended over circular microfabricated holes. 

Raman measurements on suspended nanostructures can be used to determine their ther- 
mal conductivity. This method has been employed to measure the thermal conductivity 
of single layer graphene''"^''^^'^^ (Fig. 13). In one experiment, a single layer of graphene is 
mechanically placed across microfabricated Si02 trenches, to remove any interaction of the 
graphene layer with the substrate. The laser source used for the Raman measurement is also 
used as the local heating probe. By monitoring the shift in the G peak as a function of the 
change in laser power P, the thermal conductivity K of a graphene layer can be obtained 
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Single layer 
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Bilayer 


2D 
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TABLE 111. Temperature coefficients for the G and 2D peaks in suspended graphene layers. The 
values of w^ are extrapolated by fitting-'^ ''^^. The x values for the G peak are compared to those 
for other carbon-based materials. 



according to 



53. 



K 



XG{L/2hW) 
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where xg is the temperature coefficient of the G peak, L is the the distance from the middle 
of the suspended graphene layer to the heat sink, h and W are the thickness and width 
of the graphene layer, respectively. Equation 12 is valid under the assumption that the 
front wave is non spherical, as is usually the case when the laser spot size (~ 0.5 — 1.0 /xm) 
is of the same order as the graphene strip lateral size'^'^ Although the interaction with 
the substrate is minimized across the trenches, residual strain may still be present in the 
supporting regions. The amount of strain in the suspended region however was considered 
negligible, as the Raman peak position in this region, at room temperature, corresponds 
to that of unstrained suspended graphene (Fig. 6) '■'''■^'*. Furthermore, the coefficient xg in 
eq. 12 is measured on an unsuspended graphene layer, while the experiment is performed on 
a suspended layer, xg foi' an unsuspended graphene monolayer is expected to be lower than 
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FIG. 13. The thermal conductivity measurement is performed by monitoring the change in position 
of the G peak as a function of the total dissipated power. The excitation laser light focused on 
a single layer graphene suspended across a trench (inset), is to create a local radiative hot spot, 
and to generate a heat wave across the graphene layer ''. (Reprinted with permission from ref.'"''^. 
Copyright 2008 American Chemical Society.) 

that for the suspended layer because of the interaction between the graphene layer and the 
substrate. Therefore, the measured thermal conductivity is underestimated. As previously 
noted, Cai et al. performed a similar experiment where graphene synthesized via chemical 
vapor deposition was pressed against a Au/Si02 thin film on Si'"'^ . The significant difference 
in the value of xg measured in these experiments may be due to an enhanced interaction of 
graphene with the substrate, possibly due to the graphene synthesis and deposition method. 
Under these conditions, the coefficient of thermal expansion of graphene is strongly affected 
by that of the substrate, leading to a value of xg which is significantly different from that 
of a purely suspended graphene film. Further investigations are needed to quantify how 
the thermal evolution of the graphene Raman spectra is affected by the graphene-substrate 
interaction and in particular by the difference in the coefficients of thermal expansions of 
graphene and the substrate. 

The measured thermal conductivity is compared to those of other carbon based materials 
in table IV. The extremely high value of phonon thermal conductivity in the strictly two- 
dimensional graphene layer is in sharp contrast with the reduced phonon thermal conductiv- 
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TABLE IV. The thermal conductivity of a single and multilayer layer graphene is measured opti- 
cally via Raman spectroscopy. Thermal conductivities of single layer graphene, single- and multi- 
wall CNT and diamond (via the 3-a; method) are showed for comparison. 



ity (as compared to bulk values) in quasi-one- dimensional systems such as nanowires""", or 
quasi-two-dimensional semiconducting thin films'^ '. The net reduction in the phonon thermal 
conductivity observed in these systems is explained in terms of rough boundary scattering 
or phonon spatial confinement effects. Given the high values measured for a single layer 
graphene, such effects appear not to be present. Furthermore, when comparing the thermal 
conductivity of a single layer graphene to other graphitic materials such as CNT, graphene 
exhibits a higher value, possibly due to a reduced number of structural defects, and a reduced 
intralayer scattering. In a comparison with bulk graphite, thermal conductivity approaches 
that of bulk as the number of atomic planes in graphene films increases from 2 to 4"". It 
has been shown that Umklapp-limited thermal conductivity of graphene grows with the in- 
creasing linear dimensions of graphene flakes and can exceed that of the basal planes of bulk 
graphite when the flake size is on the order of a few micrometers^^ '^^. 
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VI. CONCLUSIONS AND OUTLOOK 

Raman spectroscopy is currently used as a metrology tool to determine the extent, the 
quality and the uniformity of graphene films. This review has illustrated the applications 
of Raman spectroscopy to probing the mechanical properties of graphene films. The direct 
measurement of Raman peak shifts, for example, has enabled the determination of param- 
eters such as the Griineisen parameter and the shear deformation potential, and thus to a 
measurement of the strain in graphene films. While such shifts, in general, can be attributed 
to other causes (e.g. induced charge, doping), under precise experimental conditions (ther- 
mal equilibrium, constant pressure, and with fixed Fermi level) lattice strain can be directly 
measured from peak changes in the Raman spectra'*''. Understanding the evolution of strain 
in graphene films is important, as it allows for a deeper understanding of how graphene 
interacts with the environment, and particularly with a substrate. The ability to monitor 
and control strain in graphene could be crucial during device fabrication, as it affects the 
electronic properties of the material itself'"*. For example, it has been recently shown that 
modulation in electrical*^^ and optical^"*^ conductance can be induced by strain. It has been 
suggested that by properly modulating strain locally in graphene may lead to a controlled 
tuning of the electronic band gap^"^. Such studies are in their infancy, however. The vast 
majority of investigations have been performed either on exfoliated graphene, or on epitaxial 
graphene grown on SiC More investigations are needed to understand the presence and the 
evolution of strain in graphene grown, for example on transition metals via chemical vapor 
deposition, or as an effect of the mechanical transfer in the case of exfoliated graphene. Since 
deposition or synthesis methods strongly affect the graphene interaction with the substrate, 
further studies are needed to highlight and establish a connection between the strength 
of this interaction and the thermal evolution of the Raman spectra of graphene. While 
attempts to correlate strain to other structural properties of graphene (such as surface mor- 
phology) have been proposed, more work is needed to be able to connect strain with the 
electrical, optical and thermal properties of the material. As doping strongly affects strain 
in thin films'"', more investigations are required to determine how doping affects the strain 
in graphene films. 

From a fundamental standpoint, Raman spectroscopy can provide accurate in situ mea- 
surements of thermal properties such as the thermal conductivity. Such approach allows for 
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the characterization of the role of geometry, chemistry, and morphology, and of their effects 
on thermal properties. Such capabilities need to be extended to other graphene-related 
materials, such as graphene oxide^°'^'^°' and graphane^^^. When applied to graphene in a 
controlled environment, these measurements, may prove suitable for sensing applications. 
Overall, the characterization of mechanical properties of graphene with Raman spectroscopy 
will promise to be valuable in the determination of the optimal growth conditions, and even 
more in the optimization of fabrication methods of graphene-based devices. 
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